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ABSTRACT 


Recordings  on  photoqraphic  Film  using  binary  data  and  usinq 
256  qrey  levels  were  made  to  simulate  signals  received  from  the 
modified  APS-9A  radar.  Included  in  the  calculation  of  the 
simulated  signal  was  the  hiqh  spatial-frequency  fall-  off 
inherent  in  film  recording. 


RESUME 


II  est  traitP  de  simulations  par  donnPes  binaires  et  donnPes 
b 256  tons  de  qris  d ' inter f Programmes  du  radar  APS-94  modofiP. 
On  a tenu  compte  de  l’amortissement  dO  h 1 * intPgrat ion  produite 
au  moment  de  l'Pcriture  sur  film. 
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INTRODUCTION 


For  some  time  there  has  been  interest  in  acquisition  and 
analysis  of  data  from  synthetic  aperture  radar  (SAR)  at  the 
Defence  Research  Establishment  Ottawa  (DREO)  (1).  An  optical 
correlator  has  been  designed  and  constructed  (2)  by  Applied 
Physics  Specialties  Ltd.  of  Toronto,  and  set  up  at  DREO.  This 
correlator  was  intended  for  producing  images  from  inter f erograms 
produced  by  a Canadian  Forces'  APS-94D  radar.  This  Argus-borne 
radar  is  being  modified  by  workers  at  the  Communications  Research 
Centre  (CRC)  to  operate  in  a SAR  mode. 

Verification  of  the  performance  of  the  correlator  requires  an 
interf erogram  of  known  and  easily  measurable  imagery.  A line  of 
points  spaced  at  predetermined  distances  would  be  an  ideal 
object.  Such  an  object  would  permit  measurement  of  the  point 
size  and  the  minimum  resolvable  distance  between  points. 
However,  considerable  effort  is  required  to  construct  and  arrange 
such  targets  on  the  ground.  Furthermore,  numerous  trials  would 
be  necessary.  The  small  amount  of  flying  time  available  to  the 
CRC  SAR  and  lack  of  staff  made  that  approach  impractical. 
Furthermore,  it  was  evident,  from  initial  difficulties  in 
converting  the  APS-94D  incoherent  radar  to  a SAR,  that  the 
optical  correlator  would  be  ready  before  the  CRC  radar.  It  was, 
therefore,  decided  to  produce  simulated  interferograms.  The 
parameters  used  for  the  simulation  were  chosen  so  as  to  describe 
as  closely  as  possible  the  characteristics  of  the  CRC  radar. 
Note  that  a simulation  has  the  advantage  of  eliminating 
electronic  noise. 


CHARACTERISTICS  OF  THE  RADAR 


The  CRC  radar  is  a modification  of  the  APS-94D  incoherent 
radar  made  by  Motorola.  The  signal  from  the  SAR  modulates  a CRT. 
The  pattern  on  the  CRT  is  recorded  on  film  which,  upon 
development,  is  the  desired  interf erogram.  The  main  parameters 
(3)  of  the  modified  recording  system  are: 


film  speed  : 1.60  mm/sec 

azimuthal  width  of  the  scanning 

beam  at  the  recording  film  plane  t 16.3  microns 

scan  rate  of  CRT  : prf  : 300  acans/sec. 

range  demagnification  factor  q : 230000 

azimuth  demagnification  factor  p : 62300 
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aspect  ratio  : K=q/p=4 

offset  frequency  of  spatial  carrier  : IS. 6 lines/mm 
bandwidth  of  signal  : 31.23  lines/mm 
wavelength  of  the  radar  : 3.3  cm 
azimuthal  mean  focal  length  : 86.8  mm 

azimuthal  angular  width  of  the  radar  beam  pattern  : .008  rad 
length  of  the  radar  pulses  : .1  micro-second 


The  radar  has  two  modes  of  operation.  First  a short  range 
mode  where  the  range  varies  from  3.3  km  to  23  km  and  a long  range 
mode  where  the  range  varies  from  23  km  to  30  km.  The  altitude  of 
flight  of  the  aircraft  is  2.28  km  (see  Fig.  1). 


A scanning  rate  of  300/sec  with  a writing  beamwidth  of  16.3 
microns  and  a film  speed  of  1.6  mm/sec  causes  a superposition 
whose  effects  will  be  discussed  later. 


SIGNAL  FROM  A POINT 


The  signal  (4)  from  a point  on  the  qround  located  at 
coordinates  (x,r)  is  described  by  i 


(xf,rf)  - b + a rect 


rf  - r/q 

a rect 

Ar/q 

cos 

2*f  xf  ♦ (*f  - jj) 


(1) 
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= azimuthal  coordinate  on  the  ground 
r s range  coordinate  on  the  ground 

Xf  s azimuthal  coordinate  on  the  film 
Tj.  = range  coordinate  on  the  film 

f s offset  freguency  of  spatial  carrier  in  lines/mm 
X s radar  wavelength 

Ar  = spatial  width  of  transmitted  pulse 

b : bias  level 

O'  * reflectivity  of  a point  reflector  at  (x,r) 


The  bias  level  has  to  be  sufficiently  large  that  the  signal 
never  becomes  negative.  The  rectangle  function  describes  the 
shape  of  a short  impulse  in  the  range  coordinate.  The  width  of 
this  pulse  sets  the  range  resolution  of  the  radar.  The  first 
term  in  the  cos  argument  is  the  azimuthal  freguency  carrier  of 
the  signal.  The  second  term  is  a quadratic  phase  factor  with  a 
range  dependent  curvature  radius  f(r)  so  that  the  azimuthal 
object  will  be  at  a distance  f(r)  from  the  film  plane.  Xe is  the 
wavelength  of  the  light  in  the  optical  correlator. 

f (r)  * \r/2p2A0  (2) 


For  construction  of  the  interf erogram  it  is  required  to 
reproduce  on  film  the  function  (1).  For  the  simulated 
inter f erograms  the  required  optical  density  variation  was 
represented  both  in  binary  form  and  with  grey  levels.  The 
simulations  with  grey  levels  have  236  levels  written  by  the 
OPTRONICS  microdensitometer  at  the  Canadian  Center  for  Remote 
Sensing  (CCRS). 


SUPERPOSITIONS  IN  THE  WRITING  PROCESS 


There  is  a smoothing  phenomenon  in  the  writing  process  due  to 
the  fact  that  the  width  of  the  writing  beam  is  16.3  microns  and 
the  speed  of  the  film  is  1.66  mm/sec.  A scanning  rate  of  100/aec 
would  just  give  16.3  microns  juxtaposed  strips.  But  the  aweepinq 
rate  is  300/sec.  Therefore,  the  film  does  not  heve  the  time  to 
advence  enough  between  two  successive  scans  to  avoid 
superposition  and  each  3.3  micron  strip  is  a superposition  of 
five  different  samples  (see  Fig.  2). 
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Successive  Positions  of 
the  16*5 Beam 


Fig.  2.  Scanning  superpositions 

The  effects  of  the  superpositions  Mere  calculated  by  means  of 
an  APL  program  by  varying  the  amount  of  superposition  in  the 
signal.  That  program  divided  each  16.5  micron  strip  in  n 
divisions  of  width  of  16.5/m  and  calculated  the  intensity  of  each 
division  by  averaging  the  intensity  of  the  n preceding  divisions. 
When  there  is  no  superposition,  that  is,  when  the  signal  is 
sampled  each  0.01  sec  and  the  resulting  16.5  micron  strips  are 
juxtaposed,  the  resulting  signal  looks  undersampled  compared  to 
the  original  signal.  When  the  number  of  superpositions 
increases,  the  undersampling  disappears  and  is  replaced  by  a 
damping  of  the  high  frequencies  (see  Figs.  2,  3,  4,  5,  6,  7 and 
8).  The  range  for  the  signals  shown  in  the  following  figures  is 
13  km. 
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Fig.  8.  Signal  with  fifteen  superpositions 


We  notice  that  the  damping  saturates  quickly  with  the 
increasing  number  of  superpositions.  There  is  almost  no 
difference  between  signal  built  from  five  and  fifteen 
superpositions.  Since  the  damping  i3  more  pronounced  for  high 
frequencies,  it  is  interesting  to  study  the  frequency  content  of 
the  signal  as  the  range  increases.  In  the  following  figures,  we 
used  five  superpositions  for  siqnals  located  at  3.3,  13,  25  and 
30  km. 
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Fiq.  12.  Signal  from  a point  object. at  a range  of  50  km 


Me  note  an  increasing  dampinq  of  the  hiqh  frequencies  as  the 
range  decreases. 

POSITION  AND  AZIMUTHAL  WIDTH  OF  TRACES 

The  azimuth  width  of  the  SAR  beam  pattern  on  the  ground  and 
the  time  position  of  the  returned  siqnal  are  both  functions  of 
the  range  which  varies  between  3.3  and  50  km.  The  angular  width 
in  azimuth  of  the  beam  is  0.008  radian  so  the  width  on  the  qround 
is  given  by 

l * 0.008  x r (3) 
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A point  on  the  ground  will  return  a signal  only  if  it  is  in 
the  beam  of  the  radar.  If  we  set  t = 0 when  the  point  on  the 
ground  passes  the  aircraft,  the  return  signal  exists  between  time 

-k  < t < fc  <*> 

where  v is  the  aircraft  speed  (100  m/sec). 

If  the  speed  of  the  film  is  1.66  mm/sec  and  if  the  position 
x=0  on  the  film  corresponds  to  the  position  t=0  of  the  aircraft, 
the  time  limits  give  the  following  position  limits. 

1.66  x l _ „ _ 1.66  x i ro 

2v  < Xf  2v 

Figures  13,  14  and  13  illustrate  some  examples  of  the  siqnal 
to  be  put  on  film. 


+ l'65mm 


Fig.  14.  Signal  for  a range  of  13030  m 


+ l-65mm 


Fig.  13.  Signal  for  a range  of  16060  m 
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Fiq.  16.  Signal  for  a range  of  24940  m 


BINARY  SIMULATIONS 


Binary  simulations  of  the  signal  we  first  constructed  to  get 
a quick  and  inexpensive  verification  for  the  correlator.  To  do 
so,  we  determined  the  location  of  the  signal's  zero  crossings  and 
had  the  draftsmen  draw,  at  100X  scale,  a series  of  black  and 
white  strips  between  the  zero  crossings.  The  drawing  was 
photo-reduced  by  100  (see  Fig.  17)  to  produce  a transparency  such 
that  a transparent  strip  represents  any  signal  qreater  than  the 
bias  level  and  an  opaque  strip  represents  signals  smaller  than 
the  bias. 


Fig.  17.  Binary  simulation  for  a range  of  19030  m 


The  simulated  object  was  a line  at  a given  range.  We  made 
simulations  for  five  different  values  of  the  renge:  13030, 
16060,  19030,  22060  and  24940  m.  Note  that  the  binary 
simulations  do  not  take  into  account  the  damping  associated  with 
the  writing  process. 
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SIMULATIONS  WITH  256  GREY  LEVELS 


We  made  a more  8ophisticated  simulation  taking  into  account 
the  damping  of  the  08cillations  and  their  continuous  character. 
We  used  the  CCRS  OPTRONICS  microdensitometer.  It  can  write 
linearly  256  grey  levels  with  a 25  micron  resolution.  The  data 
must  be  given  to  the  microdensitometer  as  density  quantified  into 
256  grey  levels.  Density  D is  defined  as 


where 

I 

I 


i 

t 


D 


log 


incident  intensity 
transmitted  intensity 


(6) 


§5 


The  microdensitometer  uses  2474  film  from  Kodak.  That  film, 
when  processed  in  good  conditions,  has  a Y =1  up  to  a density  of 
2.3.  So  we  must  add  to  the  original  signal  vsrying  between  -1 
and  +1  the  bias  level  b and  a scale  factor  to  obtain  density 
varying  between  0 and  2.3.  After  that,  we  have  to  take  the 
logarithm  and  to  quantify  the  densities  in  256  levels  (see  Fig. 
18). 


+ l-65mm 


Fig.  18.  Density  distribution  for  a range  of  13000  m 


The  simulated  object  is  a line  of  18  dots  spread  in  ranqe  on 
40  positions  and  spaced  so  we  can  evaluate  the  system  resolution. 
The  pattern  of  40  positions  is  repeated  periodically  with 
increasing  values  of  the  range.  Figure  19  shows  the  arrangement 
of  the  18  dots  in  the  40  positions. 
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15  10  15  20  25  30  35  40 


Fig.  19.  Arrangement  of  object  points 


The  azimuthal  values  of  intensities  are  calculated  each  3.3 
microns  and  written  by  the  microdensitometer  with  a 23  microns 
resolution.  A length  of  23  microns  on  the  film  produced  by  the 
microdensitometer  corresponds  to  3.3  microns  on  the  film 
written  by  the  radar.  So  it  must  be  photo-reduced  of  7.6  . The 
resolution  cell  range  width  d on  the  film  written  by  the  radar 
is 


(7) 


where 

Ar  a spatial  width  of  the  range  pulse  (120  microns) 
c = speed  of  light 
T = length  of  a pulse 


Thus  the  width  resolution  cell  in  range  on  the  film  d 
written  by  the  microdenaitometer  must  bet 


dm  * df  X 7.6  (8) 


That  is  approximately  33  times  23  microns.  So  we  must  repest 
the  information  33  times  in  range  to  get  the  right  scaling  of  the 
azimuth  relative  to  the  range.  The  film  printed  and  processed 
must  then  be  reduced  by  7.6  (Fig.  20). 
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CONCLUSION 


We  here  discussed  the  characteristics  oF  the  siqnals 
generated  by  the  modified  APS-94  radar  and  described  the 
techniques  used  to  simulate  those  siqnals.  We  made  binary 
simulations  for  5 different  range  values  and  a qrey  level 
simulation  Mhere  the  range  varied  between  13030  m and  19030  m. 
The  results  of  testa  made  on  the  correlator  using  these 
simulations  will  be  given  in  a subsequent  report  alonq  with  the 
results  from  other  testing  methods. 


UNCLASSIFIED 


UNCLASSIFIED 


IS 


REFERENCES 


1.  R.T.  Lowry,  DREO  Technical  Note  No.  74-6,  A Feasibility  Study 
for  Real  Time  Optical  Correlator  for  the  APS-94D  (U). 

2.  Report  on  Advanced  Optical  Desiqn  of  Radar  Correlator, 
Canadian  Instrumentation  and  Research  Ltd.  April  1975,  D.S.S. 
contract  No.  2SR4-0123. 

3.  E.B.  Felstead,  Personnel  Communication. 

4.  W.H.  Lee  and  M.Q.  Greer,  Matched  Filter  Optical  Processor 
Appl . Opt  4 vol.13  p 925. 


UNCLASSIFIED 


llnclaaalflpri — 

Security  Classification 


DOCUMENT  CONTROL  DATA  - R ft  0 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  document  is  classified) 

1.  originating  activity 

Defence  Research  Establishment  Ottawa 
National  Defence  Headquarters 
Ottawa.  Ontario KlA  0Z4 

2a.  DOCuy^NY^ECU^T^gLASSlF  ICATION 

2b.  GROUP 

HA 

3.  DOCUMENT  TITLE 


Simulated  Inter ferograms  of  the  Modified  APS-94D  Motorola  Radar  (U) 

4.  DESCRIPTIVE  NOTES  (Type  of  report  end  inclusive  dates) 


5.  AUTHOR(S)  (Lest  name,  fust  name,  middle  initial) 


N.  Brousseau 

6.  DOC"*-p*'t  DATE 

7a  TOTAL  OP  PAGES 

7b.  NO.  O^  REFS 

DECEMBER  1976 

8a.  PROJECT  OR  GRANT  NO. 

9a  ORIGINATOR'S  DOCUMENT  NUM8ERISI 

38-80-07 

DREO  TN  - 76-32  ^ 

8b.  CONTRACT  NO. 

9b.  OTHER  DOCUMENT  NO.(S)  (Any  other  numbers  that  may  be 
assigned  this  document) 

10.  DISTRIBUTION  STATEMENT 

Unlimited  - See  attached  list 

11  SUPPLEMENTARY  NOTES 

12.  SPONSORING  ACTIVITY 

DREO 

i3.  abstract  _ Unclassified 

Recordings  on  the  photographic  film  using  binary  data  and  using 
256  grey  levels  were  made  to  simulate  signals  received  from  the  modified 
APS-94D  radar.  Included  in  the  calculation  of  the  simulated  signal  was 
the  high  spatial- frequency  fall-off  inherent  in  film  recording. 


KEY  WORDS 


- synthetic  aperture  radar 

- simulation  of  radar  data 

- optical  processing 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  th„  name  and  address  of  the 
organization  issuing  tha  document. 

2a.  OOCUMENT  SECURITY  CLASSIFICATION  Enttr  tha  overall 
security  classification  of  tha  document  including  special  warning 
terms  whenever  applicable. 

2b.  GROUP:  Enter  security  reclassification  group  number.  The  three 
groups  are  defined  in  Appendix  'M'  of  the  ORB  Security  Regulations. 

3.  OOCUMENT  TITLE:  Enter  the  complete  document  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified.  If  a 
sufficiently  descriptive  title  cannot  be  selected  without  classifi- 
cation, show  title  classification  with  the  usual  one-capital-letter 
abbreviation  in  parentheses  immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  Enter  the  category  of  document,  e.g. 
technical  report,  technical  note  or  technical  letter.  If  appropri- 
ate. enter  the  type  of  document,  e.g.  interim,  progress, 
summery,  annual  or  final.  Give  the  inclusive  dates  when  a 
specific  reporting  period  is  covered. 

5.  AUTHORISE  Enter  the  name(s)  of  authorls)  as  Aown  on  or 
m the  document.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank.  The  name  of  the  principal  author  is  an 
absolute  minimum  requirement. 

6.  DOCUMENT  DATE:  Enter  the  date  (month,  year)  of 
Establishment  approval  for  publication  of  the  document. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count  should 
follow  normal  pagination  procedures,  i.e.,  enter  the  number 
of  pages  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  document. 

SB.  PROJECT  OR  GRANT  NUMBER:  If  appropriate,  enter  the 
applicable  research  end  development  project  or  grant  number 
under  which  the  document  was  written. 

8b.  CONTRACT  NUMBER:  If  appropriate,  enter  the  applicable 
number  under  which  the  document  was  written. 

9e.  ORIGINATOR'S  DOCUMENT  NUMBER (S):  Enter  the 
official  document  number  by  which  the  document  will  be 
identified  and  controlled  by  the  originating  activity.  This 
number  must  be  unique  to  this  document. 


9b.  OTHER  DOCUMENT  NUMBER(S):  If  the  document  has  been 
assigned  any  other  document  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s). 

10.  DISTRIBUTION  STATEMENT:  Enter  any  limitations  on 
further  dissemination  of  the  document,  other  then  those  imposed 
by  security  classification,  using  standard  statements  such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
document  from  their  defence  documentation  center." 

(2)  "Announcement  and  dissemination  of  this  document 
is  not  authorized  without  prior  approval  from 
originating  activity." 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explanatory 
notes. 

12.  SPONSORING  ACTIVITY:  Enter  the  name  of  the  departmental 
project  office  or  laboratory  sponsoring  the  research  and 
development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a brief  end  factual 
summary  of  the  document,  even  though  it  may  also  appear 
elsewhere  m the  body  of  the  document  itself.  It  is  highly 
desirable  that  the  abstract  of  classified  documents  be  unclassi- 
fied. Each  paragraph  of  the  abstract  shall  end  with  an 
indication  of  the  security  classification  of  the  information 

in  the  paragraph  (unless  the  document  itself  is  unclassified) 
represented  as  ITS).  (S).  (Cl.  (R).  or  (U). 

The  length  of  the  abstract  should  be  limited  to  20  single-spaced 
standard  typewritten  lines;  7*4  inches  long. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms  or 
short  phrases  that  characterize  a document  end  could  be  helpful 
in  cataloging  the  document.  Key  words  should  be  selected  so 
that  no  security  classification  is  required.  Identifiers,  such  as 
equipment  model  designation,  trade  name,  military  project  code 
name,  geographic  location,  may  be  used  as  key  words  but  will 
be  follow  ^d  by  an  indication  of  technical  context. 


BUREAU  - RECHERCHE  ET  DEVE LOPPEMENT 

MINISTERE  DE  LA  DEFENSE  NATIONALE 
CANADA 

CENTRE  BE  RECHERCHE*  POUR  U DEFENSE,  OTTAWA 

NOTE  TECHNIQUE:  N°  76  32 


• » 

I 


SIMULATIONS  D'INTERFEROGRAMMES  DU  RADAR 
MOTOROLA  APS-94  MODIFIE 

par 

N.  Brousseau 

Section  de  teledetection 
Division  des  sciences  de  la  terre 


A 


REgU  DECEMBRE  1976 
PUBLIE  JANVIER  1977 
OTTAWA 


PROJECT  N° 
3890-07 


11-49% 


precedhc  page^blank-not  fumed 

V 

NON-CLASSIFIF 


RESUME 


II  est  traitd  de  simulations  par  donndes  hinaires  et  donndes 
5 256  tons  de  oris  d ' interf droqrammes  du  radar  APS-94  nodifid. 
On  a tenu  compte  de  1 ' amort issement  dO  & 1 ' intdnr at  ion  produite 
au  moment  de  l'dcriture  sur  film. 


ABSTRACT 


Recordinqs  on  photoqraphic  film  usinn  binary  data  and  usinn 
256  qrey  levels  were  made  to  simulate  siqnals  received  fro*  the 
modified  APS-94  radar.  Included  in  the  calculation  of  the 
simulated  siqnal  was  the  hiqh  spat i a 1 -f r equenc y fall-off  inherent 
in  film  recordinq. 
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INTRODUCTION 


Le  Centre  de  Recherche  oour  la  Defense  d'Ottawa  (CPD0) 
s'interesse  depuis  lonqtemps  b 1 ' acquisit ion  et  b l'annlyse  (!) 
de  donnles  provenant  de  radars  & ouvertur^  synthetique  (SAR) 

Or  le  Centre  de  Recherche  sur  les  Communications  (CRC)  entrenrit 
en  1973  la  construction  d'un  SAR.  CRDO  en  a orofite  oour  mettre 
sur  pied  un  projet  parallfele  visant  & construire  le  corr61ateur 
optique  n^cessaire  au  trsitement  des  donndes  de  ce  SAR.  Ce 
corrdlateur  est  de  conception  et  de  fabrication  entiferement 
canadiennes  ( 2 ) . 


La  verification  du  bon  f onct ionnement  du  correiateur 
necessite  des  essais  avec  un  inter f er oqr amme  dont  l'imaqerie 
rdsultante  est  connue  et  facilement  mesurables.  Une  linne  de 
points  situ£8  b une  distance  predetermines  les  uns  des  autres 
serait  un  objet  ideal.  Elle  permettrait  de  mesurer  la  nrandeur 
des  points  oroduits  ainsi  que  l'espacement  minimum  permettant  de 
r6soudre  deux  points  voisins.  La  realisation  physique  d'un  tel 
objet  pose  plusieurs  oroblfemes.  L ' i nstal 1 at  ion  de  reflecteurs 
sur  le  terrain  suppose  la  disponibilite  de  la  main  d'oeuvre 
ndcessaire  h la  construction  des  cibles  ainsi  qu'b  leur  mise  en 
place.  Le  choix  du  materiau  pour  les  cibles  constitue  un  autre 
probifeme  deiicat.  En  effet,  un  materiau  possedant  un  trop  fort 
coefficient  de  reflexion  produirait  un  siqnal  pouvant  facilement 
depasser  les  limites  de  linearite  du  systfeme  d ' enr eq ist rement  et 
obliterer  une  qrande  surface  du  film.  Par  contre  un  materiau 

avec  un  coefficient  de  reflection  trop  faible  ne  serait  pas 

detecte.  La  mise  au  point  d'une  telle  experience  necessiter ait 
done  de  nombreux  essais  et  une  mobilisation  importante  de 
ressources.  Or,  la  trfes  petite  quantite  d'heures  de  vol  dont 
dispose  le  qroupe  de  CRC  ainsi  cue  la  quantite  trfes  limitee  de 
personnel  disponible  rend  cette  aporoche  impossible.  De  plus,  le 
qroupe  du  CRC  a rencontre  un  qrand  nombre  de  difficultes  dans  la 
realisation  de  son  projet  et  il  s'est  avere  assez  tfit  que  le 
corrdlateur  optique  serait  prfet  b fonctionner  avant  que  le  radar 
ne  soit  complfetement  mis  au  point.  C'est  pourquoi  nous  avons 
decide  de  remplacer  1 ' inter f eroqramme  nroduit  par  le  radar  oar 
une  simulation  de  cet  inter f eroqramme . La  simulation  possfede,  en 
outre,  l'avantaqe  d'eiiminer  la  plupart  des  sources  de  bruit 
* introduites  par  1 ' apparei 1 1 aqe  e lect ronique . II  va  de  soi  que 

les  paramfetres  de  la  simulation  ont  ete  choisis  de  fa<;on  b rendre 
le  mieux  possible  les  caracter ist iques  du  radar  du  CRC. 
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CARACTERISTIQUES  OU  RADAR  A SIMULER 


Le  radar  construit  par  le  CRC  est  une  modification  du  radar 
incoherent  APS-94  de  Motorola.  Nous  en  souliqnerons  ici  (3)  les 
principales  ca racter i st  iques : 

vitesse  de  d£roulement  du  film  : 1.60  mm/sec 

larqeur  azimutale  du  faisceau  de  balayaqe  : 16.3  microns 

frequence  de  balayaqe  : 300/sec 

facteur  d'echelle  q (nortec)  : 230000 

facteur  d'echelle  p (azimut)  : 62300 

rapport  K des  facteurs  d'echelle  : K=q/p=4 

frequence  spatiale  porteuse  : 13.6  liqnes/mm 

larqeur  de  bande  du  siqnal  : 31.23  liqnes/mm 

lonqeur  d'onde  du  radar  t 3.3  cm 

lonqeur  focale  azimutale  moyenne  : 86.8  mm 

larqeur  anqulaire  azimutale  du  faisceau  : .008  radian 

lonqeur  dea  impulsions  t .1  micro  seconde 

Ce  radar  possfede  deux  modes  de  f onctionnement : un  mode  4 
courte  portee  ou  la  portee  varie  entre  3.3  km  et  23  km  et  un  mode 
b lonque  portee  ob  la  oortee  varie  entre  23  km  et  30  km. 
(.'altitude  de  vol  de  l'avion  est  de  2.28  km.  On  a done  la 
conf iquration  suivante: 

altitude 


2*28  Km 


fiq.  1 Modes  d'operation  du  radar. 
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Notons,  en  passant,  qu'une  Frequence  de  balayaqe  de  500/sec 
avec  une  largeur  de  faisceau  d'dcriture  de  16.5  microns  et  une 
vitesse  de  ddroulement  du  Film  de  1.6  mm/sec  produit  une 
superposition.  Les  eFFets  d ' intdqrat ion  resultant  seront  dtudids 
en  ddtail  plus  tard. 


SIGNAL  A SIMULER 


Le  signal  produit  par  un  point  (4)  situd  sur  terre  aux 
coordonndes  (x  ,r  ) est  ddcrit  par  : 


ou 


Uf,ry)  = b + o rect 


rf  - r/q 
Ar/q 


cos 


(1) 


x : coordonnde  azimutale  sur  le  terrain 
r : coordonnde  de  la  portde  sur  le  terrain 
Xj.  s coordonnde  azimutale  sur  le  Film 
r^ t coordonnde  de  la  portde  sur  le  Film 
F : Frequence  spatiale  porteuse  en  liqnes/mm 
‘X  : longeur  d'onde  du  radar 

Ar:  largeur  spatiale  de  l'impulsion  transmise 
b :Fond  continu 

O'  : rdFlectivitd  d'un  point  A (x,r)  sur  le  terrain 


La  quantity  t est  un  Fond  continu  qui  doit  ttre  suFFisamment 
grand  pour  que  t(x,r)  ne  devienne  jamais  ndqatiF  4 cause  des 
oscillations  du  second  terme.  La  Fonction  rectanqle  ddcrit  la 
Forme  d'une  impulsion  courte  dans  la  coordonnde  de  la  portde.  La 
largeur  de  cette  impulsion  determine  la  resolution  en  portde  du 
radar.  Le  premier  argument  du  cosinus  est  associd  A la  porteuse 
du  signal.  Notons  que  cette  porteuse  a dt 6 choisie  selon  l'axe 
azimutal.  Le  second  terme  ddcrit  un  Facteur  de  phase  sphdrique 
d'un  rayon  de  courbure  F(r)  ddpendant  de  la  portde.  L'objet 
azimutal  sera  done  situd  A une  distance  F(r)  en  avant  du  Film.X0 
est  la  longeur  d'onde  de  la  lumiAre  utilisde  dans  le  corrdlateur 
optique. 

£(r)  - Ar/2P2X0  (?) 


Nous  nous  eFForcerons  de  reproduire,  telle  que  le  radar 
l'dcrit  aur  Film,  la  Fonction  ddcrite  par  l'dquation  1).  Pour  ce 
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faire,  nous  r£aliserons  deux  types  de  simulation.  Oans  le 
premier  cas  l'intensitd  sur  le  Film  variera  de  Fagon  binaire. 
Dans  le  second  cas  l'intensitd  sur  le  Film  comptera  256  tons  de 
gris  grflce  6 1 ' utilisation  du  microdensitomfctre  OPTRONICS  du 
Centre  Canadien  de  Tdldddtect ion  (CCT)  . 


EFFETS  DES  SUPERPOSITIONS  LORS  DE  L’ECRITURE  PAR  LE  RADAR. 


Le  phdnom&ne  d ' intdgrat ion  dont 
l'existance  provient  de  la  combinaison 
La  largeur  du  Faisceau  d'dcriture  sur  1 
et  la  vitesse  de  ddroulement  du  Film 
d£duit  que  le  taux  de  balayage  donnant 
tout  juste  juxtaposes  serait  de  100  ba 
de  balayaqe  est  de  500  /sec.  II  y a 
Film  n'a  pas  le  temps  d'avancer  suFFisa 
successiFs.  En  consequence,  chaque  tra 
une  superposition  de  cinq  intensity 
diFFdrentes  (voir  Figure  2). 


nous  avon 
su ivante 
e Film  est 
est  de  1.66 
des  bandes 
layaqes/sec 
done  super 
mment  entre 
nche  de  3. 
s voisines 


s ddjb  siqnaie 
de  paramfetres. 
de  16.5  microns 
mm/sec.  On  en 
de  16.5  microns 
. Or  le  taux 
position  car  le 
deux  balavaqes 
3 microns  est 
et  16nferement 


Positions  Successivesdu 
Faisceau  de  16*  5 

3*3/t  3*3 fi  3*3/4 


niem®  balayage 


Fig.  2 Superposition  des  balayages. 


Nous  avons  6tudi6  lea  eFFeta  de  ces  superpositions  au  moyen 
d'un  programme  dcrit  en  APL  permettant  de  verier  la  quantity  de 
superposition  du  signal  Final.  Ce  programme  calcule  l'intenaitd 
de  chaque  bands  en  Faiaant  la  moyenne  de  1' intensity  des  n 
bandes  prdeddentes.  Loraqu'il  n'y  a pas  de  superposition, 
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c ' est-4-di re  , lorsque  le  signal  est  6chant i 1 lonn6  5 tous  les  .01 
sec  et  que  les  bandes  de  16.5  microns  sont  .juxtaposes , le 
signal  resultant  aemble  sous-dchant i 1 lonn£  par  rapport  au  signal 
original.  Lorsque  le  nombre  de  superpositions  auqmente,  l'effet 
de  sous-dchantillonnage  disparalt  pour  faire  place,  tel 
qu'illustrd  par  lea  Figures  2,  3,  4,  5,  6,  7 et  8,  h un 
amort issement  des  hautes  frequences  de  plus  en  plus  prononce. 
Les  signaux  illustrds  dans  les  figures  suivantes  ont  et 6 calculus 
pour  une  valeur  de  la  portde  de  13  km. 


-•82  mm 


•82mm 


Fig.  3 Signal  original. 


+ -82mm 


Fig.  4 Signal  sans  superposition. 


Fig.  5 Signal  avec  deux  superpositions. 
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Fig.  6 Signal  avec  trois  superpositions. 


+ *82mm 


Fig.  7 Signal  avec  cinq  superpositions. 




I 

Fig.  8 Siqnal  avec  guinze  superpositions. 

On  observe  que  la  valeur  de  1 ' amort issement  semble  aaturer 
trie  rapidement  avec  le  nombre  de  superpositiona.  II  n'y  a quire 
de  difference,  en  effet,  entre  un  siqnal  construit  avec  cinq  ou 
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Fiq.  12  Signal  produit  par  un  point  situd  & une  portde  de  50km. 


On  constate  qu'un  amort issement  plus  orononcd  du  siqnal  est 
associd  I une  augmentation  de  la  frequence  du  siqnal  et  & une 
diminution  de  la  portde. 


POSITION  ET  LARGEUR  AZIMUTALE  OES  TRACES 


La  largeur  de  la  bande  dclairde  sur  le  terrain  ainsi  que  la 
position  dans  le  temps  du  siqnal  retournd  par  le  radar  sont 
toutes  deux  Fonction  de  la  portae  et  elle  peut  varier  de  3.3  & 

50km.  Puisque  le  faisceau  possfede  une  largeur  anqulaire 
azimutale  de  .008  radian,  la  largeur  dclairde  sur  le  terrain  pour 
une  portae  donnde  est  ddcrite  par  l'dquation: 

l * 0.008  x r ( 3) 

verticale 


65mm 


» 


"T,  r ■ 
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Le  point  d'azimut  0 et  de  portde  r ne  retourne  un  siqnal  nue 
s'il  est  dclaird  par  le  faisceau.  Si  l'on  pose  que  t=0  lorsque 
l'avion  croiae  le  point  rdflecteur,  le  aiqnal  de  retour  existe 
entre 


■ k < 1 < k '<*> 


ou  v est  la  vitesae  de  l'avion  (100  m/sec). 


Puiaque  le  film  se  ddroule  & la  vitesse  de  1.66  mm/sec  et  que 
la  position  x=0  du  film  correspond  & la  position  t = 0 de  l'avion, 
ces  limites  de  temps  correspondent  aux  positions: 

1 .66  x £ „ ..  . 1 .66  x £ ( 5 ) 

2v  *f  2v 


Lea  figures  13,  14  et  15  illustrent  quelques  exemples  du  siqnal 
4 mettre  sur  film  pour  diffdrentes  portdes. 


Fig.  15  Signal  pour  une  portde  de  16060  m 
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SIMULATONS  B INAIRES 


Nous  avons  d'abord  r6alis6  des  simulations  binaires  du  signal 
af in  d'obtenir  rapidement  et  & peu  de  frais  un  premier  instrument 
de  verification  du  corrilateur  optique.  Pour  ce  faire,  nous 
avons  calculi  les  intersections  du  siqnal  avec  la  valeur  0 et 
nous  avons  fait  exicuter  par  le  dessinateur,  avec  un 
gr andia8ement  100,  une  suite  de  bandes  paralliles 
alternat ivement  claires  et  opaques.  Le  debut  et  la  fin  de  chaque 
bande  correspondaient  & deux  intersections  successives,  Le 
dessin  ainsi  obtenu  fut  riduit  photoqraphiquement  d'un  facteur 
100  (voir  figure  17). 


Fig.  17  Simulation  binaire  pour  une  portie  de  19030  mitres. 

L'objet  dont  on  a simuli  1 ' interf iroqr amme  consistait  en  une 
ligne  de  points  aituia  & une  mime  portie.  Nous  avons  realise  des 
simulations  pour  cinq  valeurs  diffirentes  de  la  portie,  soient, 
13030,  16060,  19030,  22060  et  24940  mitres. 

Notons,  en  passant,  que  cette  methods  de  simulation  ne  tient 
pas  compte  de  1 ' amort isaement  dO  au  procidi  d'icriture  du  radar. 


4 
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SIMULATION  AVEC  256  TONS  DE  GRIS 


Nous  avons  rdalisd  une  simulation  plus  sophistiqude  tenant 
compte  de  1 ' amort issement  des  oscillations  ainsi  que  de  la  nature 
continue  des  variations  du  siqnal.  Pour  ce  faire,  nous  avons 
utilisd  le  microdensi tomfet re  OPTRONICS  du  CCT.  II  peut  dcrire 
lindairement  256  niveaux  avec  une  resolution  de  25  microns.  Les 
donndes  utilisdea  pour  la  construction  d'une  imaqe  doivent  dtre 
fournies  au  microdensitomfet re  sous  forme  de  densitd  quantifide  en 
256  niveaux  et  la  densitd  optique  est  ddfinie  commme  suit 


ou: 


D 


log 


(6) 


I^s  intensitd  lumineuse  incidente 
1^.=  intensitd  lumineuse  transmise 

Le  microdensitomfetre  dcrit  sur  du  film  2474  de  Kodak.  Ce 
film,  lorsqu'il  est  ddveloppd  dans  de  bonnes  conditions  ordsente 
un  y si  jusqu'b  une  densitd  de  2.3.  II  faut  done  ajouter  au 
siqnal  variant  entre  -1  et  +1  un  fond  continu  et  un  facteur 
d'dchelle  approprid  pour  obtenir  des  densitds  variant  entre  0 et 
2.3.  Nous  devons  ensuite  prendre  le  logarithme  et  quantifier  les 
densitds  en  256  niveaux  (voir  figure  18). 


Fig.  18  Densitd  pour  une  portde  de  13  km. 


L'objet  choisi  pour  cette  simulation  est  une  ligne  formde  de 
18  points  rdpartis  sur  40  positions  et  espacds  de  fagon  4 pouvoir 
dvaluer  la  rdsolution  du  systems.  Ce  patron  de  40  positions  se 
rdptte  pdr iodiqusmant  avec  lea  valeura  croissantes  de  la  portde  . 


NON-CLASSIFIE 


12 


NON-CLASSIFIE 


La  figure  19  illustre  la  di8position  des  18  points  dans  les  40 
positions. 


w//Mnw//mmw0m 


I 5 10 
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Fig.  19  Disposition  des  points  objet. 


Les  valeurs  de  l'intensite  calculdes,  en  azimut,  A tous  les 
3.3  microns,  seront  dcrites  avec  le  microdensitomfetre  fe  une 
resolution  de  23  microns.  Une  longeur  de  23  microns  sur  le  film 
produit  par  le  microdensitomfetre  correspond  done  fe  une  lonqeur  de 
3.3  microns  sur  le  film  produit  par  le  radar.  II  faudra  done 
rdduire  d'un  facteur  7.6. 


La  largeur,  selon  la  portee,  de  la  cellule  de  resolution  sur 
le  film  ecrit  par  le  radar  est  de 

d . £1  = cT  (7) 

f q g 


ou 

Ar:largeur  spatiale  de  l'impulsion  selon  la  portee 
civitesse  de  la  lumifere 
Tsdurde  d'une  impulsion 

La  largeur  selon  la  portee  de  la  cellule  de  resolution  sur  le 
film  ecrit  par  le  microdensitomfetre  devra  fttre  dqale  fet 

dm  " df  X 7,6  *n) 

Ce  qui  est  approximativement  egal  fe  33  foia  23  microns.  II 
faudra  done  rdp6ter  1 ' inf ormat ion  azimutale  33  foia  pour  obtenir 
le  bon  facteur  d'dchelle  entre  l'azimut  et  la  portee.  Une  foia 
le  film  produit,  une  photoreduction  d'un  facteur  7.6  sera 
effectude  (voir  figure  20). 


NON-CLASSIFIE 


14 


NON-CLASSIF IE 


CONCLUSION 


Nous  avons  discute  les  caractdr istiques  des  siqnaux  q6n6r£s 
par  le  radar  APS-94  modifie  et  nous  avons  d£crit  les  techniques 
utilis£es  pour  simuler  ces  signaux.  Nous  avons  r£alis£  de  la 
sorte  des  simulations  binaires  pour  diffdrentes  valeurs  de  la 
portae  et  une  simulation  en  tons  de  gris  ou  la  portde  variait  de 
13030  m & 19030  m.  Les  rdsultats  des  tests  effectuds  & l'aide  de 
ces  simulations  seront  inclus  dans  un  raoport  d^crivant  les 
diffdrentes  verifications  auxguelles  le  corrdlateur  a 6 td  soumis. 
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